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This paper describes the Computational Fluid Dynamics (CFD) model developed to 
simulate the supersonic rocket exhaust in an entrained flow cylinder. The model can be used 
to study the plume-induced environment due to static firing tests of the Taurus-II launch 
vehicle. The finite-rate chemistry is used to model the combustion process involving rocket 
propellant (RP-1) and liquid oxidizer (LOX). A similar chemical reacting model is also used 
to simulate the mixing of rocket plume and ambient air. The model provides detailed 
information on the gas concentration and other flow parameters within the enclosed region, 
thus allowing different operating scenarios to be examined in an efficient manner. It is 
shown that the real gas influence is significant and yields better agreement with the theory. 
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Nomenclature 

pre-exponential factor for reaction step k 
molar concentration of species / 
coefficient of diffiisivity 
total energy per unit volume 

activation energy of reaction step k 
fluxes in x, 3 /, z directions 
ratio of specific heats 
Gibbs free energy of species / 
coefficient of thermal conductivity 
forward rate constant for reaction k 

backward rate constant for reaction k 
first coefficient of viscosity 
second coefficient of viscosity 
chemical symbol of species i 

exponent of temperature in the rate constant of reaction k 

dependent conservation variable 

gas constant 

universal gas constant 

turbulence transport quantities 

density 
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source term 
static temperature 
viscous stresses 

velocities in x,y,z directions 
molecular weight of i 

stoichiometric coefficient of / in reactant side of step k 
stoichiometric coefficient of i in product side of step k 
mass flow rate per unit volume of production of i from reaction k 


I. Introduction 

N ASA awarded to Orbital Sciences Corperation (OSC) a Commercial Orbital Transportation Service (COTS) 
contract to demonstrate delivery of cargo to the International Space Station. For these COTS missions OSC 
intends to use Taurus-II to launch its Cygnus spacecraft. The baseline Taurus-II is a two-stage vehicle, designed to 
launch payloads weighing up to 15,000 lbs into low-Earth orbit (Figure 1). The first stage uses RP-1 (kerosene) and 
liquid oxygen (LOX) as propellants, powering two Aerojet AJ-26-62 main engines. The second stage is a CASTOR 
30 solid rocket motor, built by Alliant Tech Systems, Inc. (ATK). This paper focuses on the first stage engine only. 
This engine is previously known as the Russian NK-33 liquid rocket engine, which provides high performance with 
light weight design, high chamber pressure, stable combustion, and a LOX rich prebumer^ Currently OSC is 
planning to conduct a series of cold and hot first stage tests at Wallops Flight Facility (WFF). During these tests the 
two main rocket engines will be enclosed in an extended cylinder that resembles the hold-down bay. As shown in 
Figure 2, there are several vents around the cylinder to allow air entrainment during the hot firing test. The 
objective of this analysis is to solve the flowfield within the enclosed cylinder and to quantify the amount of air 
entrained during the test. The solutions can also be applied toward the study of the plume-induced environment. 
The analysis is based on multi-species, chemical reacting simulations. The results at the nozzle exit are compared to 
theoretical rocket performance data obtained from chemical equilibrium with applications (CEA)^. A simplified, 
single-species simulation is also presented here for comparison purposes 

II. Modeling Approach 

The analysis pertains to a complex system involving swirling flows, complex chemical reactions and complex 
geometry. These complexities provide a significant computational challenge. All modeling reported in this paper 
was performed using CFD++ 10.1, a multi-purpose CFD code developed by Metacomp Technologies, Inc.^’"^’^. 

A multi-block, structured mesh is created to model the AJ-26 rocket engines in an enclosed cylinder (Figure 3). 
There are a total of 321,936 cells for the single-species case and 406,656 cells for the multi-species case. In both 
cases a symmetry boundary condition is used to reduce the computational domain by half The entire domain is 
initialized with a small disturbance (F= 0.01 m/s); this disturbance helps initialize the freestream turbulence 
quantities in the k-s model. 

The solution is based on compressible, real gas, Reynolds-Average Navier-Stokes (RANS) with a k-s 
turbulence model. The simulation starts from the combustion chamber, where the total pressure and temperature 
are specified. For single species, ideal gas simulation only total conditions (P = 162 atm and T = 3814 K) are 
required in the combustion chamber and the gas flow in the nozzle is assumed to be air with a constant specific heat 
ratio (/=1.4). For multi-species simulations the fuel-oxidizer reactions in the chamber and mixing reactions 

outside of the nozzle are also described in the boundary conditions. The three main chemical molecules that made 
up the fuel and oxidizer are carbon, hydrogen and oxygen. Outside of the nozzle, nitrogen is added to make up the 
ambient air. Combustion is handled using 11 -species, 18-reaction chemistry mechanism solved with a finite-rate 
kinetics scheme and a constant-pressure combustion model. Standard non-reflecting, interpolated boundary 
conditions are imposed at the outflow boundary, as well as no-slip conditions at solid surfaces and symmetry 
boundary condition at the midplane. Steady-state solutions are obtained using a multi-grid, time-stepping scheme. 
The solutions are compared between single-species, non-reacting and multi-phase, finite-rate chemical reacting 
flows. The solutions are also compared against theoretical performance data obtained from CEA. 
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III. Governing Equations 

The CFD-H- flow solver can be used to provide steady-state or unsteady flowfield solutions by solving the 
transport equations such as mass conservation equation, time averaged Navier-Stokes equations, energy equation 
and other scalar transport equations. The general form of these equations can be written as: 

gg , ^U--/v) , ^(g<-gy) I o 

dt dx dy dz 

where the subscripts i and v denote the inviscid and viscous flow terms, respectively. For the Reynolds- 
averaged Navier-Stokes equations the dependent quantities and the inviscid fluxes can be written as: 
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where ct^ ’s represent turbulence kinetic energy and “undamped” eddy viscosity in the pointwise turbulence 
models. Flow with multiple species can also be computed within this framework. The first five rows represent 
terms from the standard Euler equations with the first being the energy equation followed by the continuity and three 
momentum equations. The equation of state couples pressure to density and temperature is the perfect gas equation 
of state (p = pRT ) which can be written in terms of the conservation variables as follows: 


p = {r-\) e-^\^{puf +{pvf +{pwf'^ 


The viscous terms are defined as follows: 
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where viscous stresses Zy ’S are defined as: 
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where the Stokes theorem for gases is assumed to hold true, thus X = —/u. The temperature can be related to 


the conservation quantities via the equation of state: 
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The source terms are written as: 


5 = (o 0 gy g^ Qi ••• 

where g,’s are the body forces which can be activated if necessary, and Qj’s are the source terms such as 
production and dissipation of turbulence. 


IV. Chemistry Model 

Hydrocarbon rocket fuels such as RP-1 consist of at-least 87 identifiable hydrocarbon based compounds^’ 
Modeling the reaction rates for each of these compounds in addition to the subset of possible intermediate 
compounds would result in a very large and, for the present study, unnecessary computational burden. Therefore the 
approach taken is to utilize the widely accepted wet-CO (H 2 O catalyzed) and thermal-NO mechanisms as described 
in the following sections, which model the thermophysics associated with conversion of the intermediate state 
compounds to their final combustion products. The initial global reaction steps to these intermediate molecular 
states are model as irreversible reactions of the RP-1 hydrocarbon chains with oxygen. The initial intermediate 
states in the combustion chamber are provided by the industry standard complex chemical equilibrium code CEA^. 

In the present study, the finite-rate reactions of hydrogen, oxygen, and carbon are considered for the chemistry 
model in the combustion chamber. This model employs 8 species and 12 elementary reactions. Outside of the 
nozzle exit, nitrogen molecule is considered for the plume -air mixing reaction, which involves 3 more species and 6 
more reactions. For a general chemical reaction k\ 

i i 

where the rate of production of species i from the reaction step k can be written as follows: 

I I 

The forward rate constant for each reaction step k is given by Arrhenius kinetics: 


^ik=^i(^ik -^'ik 
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Any kinetic step, including reactions with third bodies in the product and/or reactant side can be included in the 
flow solver through inputs of Np 

The backward rate constant is computed from the equilibrium condition: 
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where the change in Gibbs free energy for reaction step k is given by: 

( 3 ) 

/=1 /=! 

A summary of chemical reactions and reaction rates used in the model is given in Table 1 . The first 12 reactions 
are related to RP-1 and LOX in the combustion chamber and the last 6 reactions describe the exhaust plume mixing 
with air outside of the rocket nozzle. These reaction equations are obtained from various sources*’ 
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V ^0^ J 
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^ Af^ ^/kmol) 

1 

H + O 2 ^ OH + 0 

1.2x10"“ 

-0.91 

69156380 

2 

H 2 + 0 ^ OH + H 

1.5x10'“ 

2.00 

31624892 

3 

0 + H 2 O ^ OH + OH 

1.5x10'" 

1.14 

72256999 

4 

OH + H 2 ^ H 2 O+ H 

1.0x10" 

1.60 

13810871 

5 

O + H + M^OH + M 

1.0x10'“ 

0.0 

0 

6 

0 + 0 + O 2 + M 

1.0x10"“ 

-1.0 

0 

7 

H + H + M^H2 + M 

9.7x10'“ 

-0.6 

0 

8 

H 2 O + M ^ H + OH + M 

1.6x10"“ 

0.0 

478375275 

9 

O 2 + H 2 ^ H + OH + OH 

7.9x10'" 

0.0 

187147514 

10 

CO + OH ^ CO 2 + H 

4.4x10“ 

1.5 

-3102421 

11 

CO + 0 + M ^ CO 2 + M 

5.3x10'“ 

0.0 

-19014967 

12 

CO + 02^ CO 2 + 0 

2.5x10'“ 

0.0 

200157545 

13 

0 + N2^^N + NO 

1.4x10'" 

0.0 

314117874 

14 

N2+02^N0 + N0 

9.4x10"" 

-2.5 

537538443 

15 

NO + 0 ^ N + NO 

1.6x10'" 

1.0 

161843947 

16 

NO + M^O + N + M 

2.3x10"“ 

-0.5 

623224476 

17 

N + OH^NO + H 

4.0x10'“ 

0 

0 

18 

C02 + N^C0 + N0 

2.0x10“' 

-0.5 

33283891 


Table 1. Reaction Rate Equations 

V. Discussion of Results 

It is shown that for the non-reacting calculations, if the K-e turbulence model is used, the swirl is predicted to 
decay too quickly and therefore recirculation zones are not captured. As a result, the single -species simulation 
shows a greater air entrainment at the ceiling, whereas in the multi -species simulation the recirculation causes flow 
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reversals, thus reducing the entrainment from the top. Figures 4 and 5 show the velocity vectors at the symmetry 
plane for single-species and multi-species, respectively. Also, particle traces in Figures 6 and 7 indicate that a lower 
entrainment rate is detected for the multi-species case, due to recirculation and flow reversals. 

Furthermore, the single-species simulation with constant specific heat ratio (/) shows an over-predicted 

solution at the nozzle exit plane. Figure 8 shows the Mach contours on the symmetry plane, where M =6.83 at the 
nozzle exit centerline; which is much higher than the CEA prediction (M =3.667). The simple explanation is that 
the speed of sound is a function of specific heat ratio and molecular weight or singularly local static temperature. 
What is expected to be occurring in reality is additional energy release from continued exothermic reaction in the 
nozzle, most notably the re-combination of dissociated molecules of hydrogen and oxygen from the now changing 
equilibrium conditions as the flow expands and reduces static pressure and temperatures in the nozzle. This is 
clearly evident in Fig 9 which depicts a sharp increase in total temperature of the flow as the flow expands past the 
throat of nozzle for the reacting chemistry case and matches well with CEA predictions. 

For single-species, ideal gas simulation, both of these parameters are assumed to be constant, resulting in a low 
speed of sound compared to real gas. Figure 10 shows the speed of sound contours for the ideal gas simulation, 
where c =384.84 m/s at the nozzle exit centerline, compared to c =879.4 m/s predicted by CEA. The under- 
prediction by the ideal gas assumption can also be found from observing other flow parameters, e.g. pressure and 
temperature at the nozzle exit centerline are predicted to be 0.0458 atm and 368.13 K, respectively. CEA yields at 
the same location a pressure of 0.664 atm and temperature of 2041.12 K. Figures 1 1 and 12 show the pressure and 
temperature contours at the symmetry plane. These results are sound in the sense that it correctly predicted the 
shock structures in the plume flowfield but the under-prediction of pressure magnitude makes the nozzle more 
over-expanded that it actually is. 

The multi-species simulation with finite-rate chemistry shows a much better agreement with the theory. Figures 
13-16 show the contours for Mach number, sonic velocity, pressure, and temperature at the symmetry plane. The 
conditions at the nozzle exit plane obtained from these simulations match up with CEA within 1% of error. In 
addition. Figures 17-21 show the mole fractions of some of the gases that made up the plume and ambient air. 
Apparently, there are strong recirculation zones between the plumes and in the region between the cylinder and 
external nozzle surface. These behaviors were not captured in the single-phase, ideal gas simulation. Table 2 
summarizes the nozzle exit plane conditions, predicted by CEA and by real gas simulation. Excellent agreement has 
been achieved with the current chemistry model. 


FLOW VARIABLES 

Chamber 

Throat 

Exit 

(CEA) 

Exit 

(CFD++) 

P, atm 

162.02 

93.541 

0.66433 

0.65199 

T,K 

3814.03 

3617.25 

2041.12 

2041.24 

Sonic velocity, m/s 

1224.2 

1182.2 

879.4 

883.7 

Mach number 

0.0 

1.0 

3.695 

3.669 

MOLE FRACTION 


CO 

0.29698 

0.28981 

0.23532 

0.23511 

C02 

0.16922 

0.18276 

0.27071 

0.27078 

H 

0.02327 

0.02030 

0.00072 

0.00072 

H2 

0.06879 

0.06675 

0.07650 

0.07689 

H20 

0.34094 

0.35551 

0.41627 

0.41570 

0 

0.01191 

0.00925 

0.00000 

0.00000 

OH 

0.06620 

0.05660 

0.00048 

0.00043 

02 

0.02238 

0.01885 

0.00001 

0.00001 


Table 2. AJ-26 Combustion Data 


VI. Conclusion 

Calculations were performed on a dual nozzle, RP-l/LOX combustion engine as part of an ongoing effort to 
study the plume-induced environment due to stage testing of Taurus-II launch vehicle. It has been shown that the 
finite-rate chemistry impacted the computational results. When included in the computational model, the chemical 
reaction significantly influences the flowfield. While the non-reacting, single-species simulation can provide 
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quantitative results, it cannot capture all physical features of this particular problem. The plume flowfield obtained 
from this model can be used to improve the thermal analysis that determines the plume impingement impact on 
ground support equipment. 
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Fig. 1 Taurus-II Configuration 
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Fig. 2 Stage 1 Test Configuration 



Fig. 3 Computational Domain 
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Fig. 4 Velocity vectors, single species 



Fig. 5 Velocity vectors, multi species 
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Fig. 8 Mach contours at the symmetry plane (single-species) 
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Fig 10 Sonic velocity contours at the symmetry plane (multi-species) 
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Fig 1 1 Pressure contours at the symmetry plane (single-species) 
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Fig 12 Temperature contours at the symmetry plane (single-species) 



Fig 13 Mach number contours at the symmetry plane (multi-species) 
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Fig 14 Sonic velocity contours at the symmetry plane (multi-species) 
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Fig 15 Pressure contours at the symmetry plane (multi -species) 
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Fig 16 Temperature contours at the symmetry plane (multi-species) 
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Fig 17 C02 mole fraction at the symmetry plane 
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Fig 18 H2 mole fraction at the symmetry plane 
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Fig 19 H20 mole fraction at the symmetry plane 
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Fig 20 02 mole fraction at the symmetry plane 
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Fig 21 N2 mole fraction at the symmetry plane 
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